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摘 要 : 以 黑河 中 游 正义 峡 水 文 站 的 实测 径流 量 数据 为 基 
检验 法 .Pettitt 检 验 法 和 累积 距 平 法 分 析 了 正义 峡 水 文 站 1970—2020 年 径流 序列 的 变化 趋势 和 突变 年 份 ,并 利用 基 
于 Budyko 假 设 的 水 热 厢 合 平衡 方程 ,对 正义 峡 径 流量 变化 趋势 进行 了 归 因 分 析 。 结 果 表 明 :(1) 研究 期 内 正义 峡 


E 础 ,采用 线性 回归 法 .Mann-Kendall 突变 检验 法 38 2/1 t 


径流 量 波动 变化 , 丰 枯 交替 ,但 整体 呈现 增加 趋势 ,2004 年 径流 发 生 突变 ,突变 后 的 年 平均 径流 量 增加 了 3.08x10* m, 
增加 率 为 32.7%。(2) 突变 后 的 2005 一 2020 年 ,正义 峡 径 流量 对 降水 . 洪 在 蒸发 和 下 垫 面 参数 的 弹性 系数 分 别 为 
1.40,-0.40 -1.57, 且 各 因子 对 径流 的 贡献 率 分 别 为 42.73% 、-12.52% .69.79% ,表明 径流 量 对 下 热 面 变化 最 为 敏感 ， 


气候 因子 中 降水 对 径流 的 影响 大 于 潜在 蒸发 .(3) 在 一 定 的 


区 域 气候 条 件 下 ,植被 覆盖 、 土 地 利用 流域 调 水 等 人 类 


活动 引起 的 中 游 下 垫 面 变化 是 正义 峡 径 流量 变化 的 主要 原 


因 。 研 究 结果 可 为 流域 管理 部 门 制定 水 资源 合理 分 配 


及 调用 方案 提供 科学 依据 。 
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在 全 球 气候 变 暖 的 背景 下 ,降水 .气温 2 
气象 因素 均 发 生 了 不 同 程度 的 变化 ,进而 导致 全 球 
许多 河川 径流 也 发 生 了 显著 变化 "“。 我 国 西 北 地 
区 对 气候 变化 的 响应 更 加 敏感 , 近 几 十 年 来 ,在 气 
温 升 高 .降水 量 增加 和 冰冻 圈 退 缩 等 众多 因素 的 影 
响 下 ,该 区 域 水 资源 时 空 分 布 更 加 不 合理 , 极 大 影 
响 了 生态 系统 稳定 。 与 此 同时 , 随 着 “西部 大 开 
发 "战略 举措 的 实施 推动 和 区 域 社会 经 济 发 展 水 平 
的 不 断 提 高 ,西北 地 区 的 土地 利用 格局 也 发 生 了 较 
大 的 改变 ,也 直接 或 间接 地 影响 着 流域 产 汇 流 过 
程 ,并 导致 一 些 河川 径流 量 的 增多 或 减少 ,这 种 变 
化 对 于 区 域 水 循环 产生 了 深远 的 影响 ”。 因 此 , 探 
究 西 北 内 陆 河 径流 变化 特征 及 影响 因素 ,对 于 深刻 
认识 干旱 区 的 生态 环境 演变 和 保障 社会 经 济 可 持 
续 发 展 具有 重要 意义 。 

目前 ,许多 学 者 对 西北 干旱 地 区 内 陆 河 流域 径 
流 变 化 及 其 影响 因素 开展 了 大 量 研究 。 有 学 者 研 
究 发 现 , 气 候 变 化 是 导致 玛 纳 斯 河和 塔里木 河 径流 
变化 的 主要 影响 因素 ,其 中 降水 增加 直接 影响 径流 
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量 增加 , 随 着 气温 增加 ,冰川 融雪 增加 ,间接 影响 玛 
纳 斯 河和 塔里木 河 径 流量 增加 “”。 也 有 学 者 对 河 
西 走廊 三 大 内 陆 河 出 山 径 流 进行 分 析 人 研究 ,结果 表 
明石 羊 河 径流 减少 ,黑河 与 踊 勤 河 径 流 呈 现 增加 趋 
势 , 三 大 内 陆 河 的 出 山 径 流 变化 受 祁连山 山区 气温 
影响 较 大 ,其 中 夏季 受降 水 影响 较 大 ,冬季 和 春季 
则 受气 温 升 高 导致 的 冰川 融雪 影响 较 大 ,山区 人 类 
活动 较 少 ,对 出 山 径 流 的 影响 较 小 “”。 通 过 对 现 
有 的 众多 研究 成 果 进 行 分 析 研 判 ,发 现 大 多 数 研 究 
主要 集中 于 内 陆 河 出 山 径流 量变 化 的 影响 因素 的 
分 析 , 这 些 人 研究 都 表明 气候 因子 和 人 类 活动 共同 对 
水 文 循环 过 程 河川 径流 演变 产生 影响 ,二 者 的 作 
用 经 常 互 相交 织 并 相互 影响 ,其 中 主导 因素 为 气候 
变化 。 内 陆 河 中 游 一 般 拥 有 适宜 的 气候 条 件 和 广 
交 的 土地 资源 ,是 工农 业 生 产 和 人 口 的 集中 分 布 
区 ,同时 也 是 水 量 消 耗 最 多 的 区 域 。 在 这 个 区 域 学 
者 们 关注 更 多 的 是 流域 调 水 、 农 业 用 水 .水 资源 承 
载 力 、. 土 地 利用 和 植被 生态 建设 “”。 中 游 也 因为 
人 类 活动 影响 频繁 ,水 文 循环 过 程 更 趋 复杂 ,径流 
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变化 研究 结果 的 时 效 性 、 可 比 性 也 欠 佳 ,对 中 游 径 
流 变化 及 其 影响 因素 的 研究 尚 显 不 足 。 因 此 ,加 强 
内 陆 河中 游 径流 演变 和 影响 因素 的 分 析 研 究 ,对 于 
揭示 该 区 域 水 文 效应 的 变化 机 理 和 如 何 应 对 水 资 
源 管 理 挑战 很 有 必要 。 

黑河 是 我 国 第 二 大 内 陆 河 ,横路 青 甘 蒙 三 省 
区 ,作为 西北 重要 的 生态 功能 区 和 生态 安全 屏障 ， 
承担 着 水 源 涵 养 .生物 多 样 性 保护 和 防风 固沙 等 诸 
多 重要 的 生态 功能 ,并 发 挥 着 极其 重要 的 生态 作 
用 。 正 义 峡 地 处 甘肃 省 张掖 市 . 酒 录 市 和 内 蒙古 额 
济 纳 旗 三 方 交界 之 地 ,是 黑河 干流 中 游 和 下 游 的 分 
界 点 ,探究 黑河 中 游 正义 峡 径 流 演变 规律 和 影响 径 
流 的 主要 因素 ,对 保障 中 下 游 生产 生活 用 水 及 维护 
区 域 生态 安全 意义 重大 。 本 文 以 黑河 中 游 长 序列 
径流 、 气 象 和 土地 利用 等 数据 资料 为 依据 ,综合 M-K 
突变 检验 法 .滑动 :检验 法 Pettitt 检验 法 .累积 距 平 
法 对 中 游 正 义 峡 径 流 演变 从 时 间 尺 度 上 进行 分 析 ， 
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并 采用 Budyko 模型 探究 径流 变化 的 主要 影响 
素 。 研 究 结果 可 为 变化 环境 下 黑河 中 游 地 区 水 资 
源 合理 开发 利用 与 科学 配置 提供 参考 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

Ri inf (98930' -101930' E, 37550' -42?40' N) 发源 
PAKÉ LI JURE , HR LETS; TALIS VA EDS EIE , 5 TA Ue 
至 正义 峡 之 间 为 中 游 ,正义 峡 以 下 为 下 游 。 黑 河中 
游 地 处 甘肃 河西 走廊 的 中 部 ,河道 全 长 185 km, DX. 
域 面 积 2.56x10' km , 属 温带 干旱 区 ,年 温差 较 大 ,年 
日 照 时 数 长 达 3000~4000 h ,年 均 温 度 7~8 *C 、 降 水 
量 110-150 mm Z$ £ Œ 1300~2400 mm .相对 湿度 
52%0 2 。 黑 河 干 流 沿 线 的 甘 州 、 临 泽 ,高台 三 县 
(区 ) 为 河西 走廊 主要 商品 粮 生产 基地 ,也 是 中 游 地 
区 河川 径流 的 主要 消耗 区 (图 1) ,近年 来 受 人 类 活 
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图 1 黑河 流域 示意 图 
Fig. 1 Map of Heihe River Basin 
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动 和 气候 变化 的 影响 ,黑河 中 游 干 流 的 年 径流 量 发 
生 了 明显 变化 。 
1.2 数据 来 源 

正义 峡 水 文 站 是 黑河 中 游 的 国家 基本 水 文 
站 。 本 文系 统 性 地 收集 了 正义 峡 水 文 站 1970— 
2020 年 的 年 径流 资料 ,来 自 甘肃 省 水 文 水 资源 局 ; 
同时 收集 了 中 游 地 区 高 台 、 临 泽 、 张 掖 三 个 气象 站 
点 1970 一 2020 年 的 气象 数据 ,来 自 国家 气象 科学 数 
据 中 心中 国 气 象 数 据 网 (http://data.cma.cn) ,包括 降 
水 湿度 .日照 时 数 ` 气温、 风速 等 资料 ; 归 一 化 植被 
指数 (NDVI) 来 自 中 国 科学 院 资 源 环境 科学 与 数据 
中 心 Chttps://www.resdc.cn) ; 潜在 蒸发 量 ET X HIK 
合 国 粮食 及 农业 组 织 (FAO) 推 荐 的 彭 曼 -蒙特 斯 公 
式 (Penman-Monteith ) 计 算 获 得 。 

1.3 研究 方法 

1.3.1 趋势 和 突变 检验 采用 线性 回归 法 .5 a 滑 动 
平均 法 对 正义 峡 水 文 站 1970 一 2020 年 径流 量变 化 
趋势 进行 分 析 ; 利用 Mann-Kendall 突变 检验 | 28 2J] t 
检验 ~ Pettitt 检验 法 和 累积 距 平 法 进行 径流 突变 检 
验 , 确 定年 径流 序列 的 突变 点 。 

(1) M-K 突 变 检 验 , 对 时 间 序 列 构建 一 个 服从 
标准 正 态 分 布 的 秩序 列 UF, 再 将 时 间 序 列 逆序 排 
列 ,对 逆序 排列 后 的 新 序列 构建 一 个 新 的 秩序 列 
UB, 则 UF 和 UB 两 条 曲线 的 交点 即 为 时 间 序 列 的 突 
变 点 。 

(2) 滑动 :检验 , 任 取 一 点 作为 基准 点 , 则 该 基 
准点 把 时 间 序 列 分 为 两 段子 序列 ,如 果 两 段子 序列 
的 均值 差异 超过 了 一 定 的 显著 性 水 平 , 则 认为 有 突 
变 发 生 , 即 该 点 为 时 间 序 列 的 突变 点 。 

(3) Pettitt 检验 ,在 连续 的 时 间 序 列 中 ,假设 一 
个 时 间 点 为 突变 点 ,检验 该 点 前 后 两 段 数据 的 累积 
分 布 函数 是 否 有 显著 差异 ,有 显著 差异 则 表明 该 点 
为 时 间 序 列 的 突变 点 。 

(4) 累积 距 平 法 ,计算 时 间 序 列 的 多 个 累积 距 

平 值 并 绘制 累计 距 平 曲线 图 ,曲线 上 升 表 示 距 平 增 
加 ,曲线 下 降 表示 距 平 减 小 ,根据 曲线 的 起 伏 变化 
可 以 判断 出 时 间 序 列 的 变化 趋势 和 发 生 突变 的 时 
间 。 具 体 计算 方法 参照 文献 [12-13] 所 述 。 
1.3.2 径流 变化 归 因 分 析 采用 基于 Budyko 假 设 
水 热 耦合 平衡 原理 的 弹性 系数 法 ”2 量化 气候 因素 
和 人 类 活动 对 正义 峡 径流 量 的 影响 。 

(1) 一 般 认为 流域 多 年 蓄 水 量变 化 可 忽略 不 


计 , 则 基于 Budyko 假 设 提出 了 流域 水 热 厢 合 平衡 方 
程 如 下 : 


pe (1) 
(P^ ET?)^ 
Px ET| 
R-P-ESP-— 2; (2) 
(P ETO" 


式 中 :及 为 径流 深 (mmy) ; P HKE (mm) ;五 为 实际 
蒸 散 发 量 (mm) ;n 为 流域 下 垫 面 参数 ; ET, 为 潜在 
ZE (mm*d") 

(2) Bi P, ET. 为 相互 独立 的 变量 ,公式 
(DAEA Rzf(P, ET, n) , 则 径流 深 的 微分 表达 
式 为 : 


-ôR oR oR 
Wt dn (3) 


(3) 用 弹性 系数 ei 表示 径流 R 对 X 因 子 的 敏感 

性 ,其 含义 为 单位 影响 因素 的 变化 引起 的 径流 变化 
量 , 其 表达 式 为 : 

gehe. (4) 


对 公式 (3) 等 式 两 边 同时 除 以 平均 径流 深 尺 并 
结合 公式 (4) 整 理 可 得 ; 
dR |... dP 
R P 
(4) 4 ec WERA AROMA 
式 (5) 可 得 降雨 (P) HERUR ET) M FRA 
件 (n) 的 弹性 系数 分 别 为 : 


dET, 。 dn 


Eo ET, "n (5) 


TE 


2 Utp e (6) 
" Qe) *6)^ - 9] 
"ne mE (7) 
+o) -( +p] 
, "In * 4^) d" Inl te”) (8) 


nil +g") -a +g) 

(5) 用 dR' 表示 基于 降雨 (P) 、 潜 在 蒸 散发 
CET, RUF Si A l'F Gn) 2E 4E, i sf Budyko [EBERT 
算出 的 突变 前 后 径流 深 的 变化 量 , dRo, dR,S. dR, 
分 别 表示 降雨 ,潜在 藻 散 发 和 下 执 面 条件 发 生变 化 
而 引起 的 径流 量变 化 量 ”: 
dR - dR, + dR, + dR, = e, AdP $ em gr ET. n a, dn 


(9) 
ICH :dP.dET, dn X P ET, 在 突变 前 后 两 段 时 间 
多 年 平均 值 的 差 值 , 则 降雨 (P) ERME ET) 
和 下 垫 面 条 件 (n) 对 径流 深 的 贡献 率 (C) 分 别 为 : 


C, - E» x100% (10) 
dR 
dR pr, 
Cpr = —À x 100% (11) 
© dR 
C, = SR. x100% (12) 
dR 
2 结果 与 分 析 


2.1 径流 变化 趋势 分 析 

2.1.1 径流 变化 趋势 ”正义 峡 水 文 站 为 黑河 中 游 干 
流 上 的 主要 控制 站 ,1970 一 2020 年 正义 峡 水 文 站 年 
径流 量 的 变化 情况 如 图 2。 由 图 2 可 以 看 出 ,正义 峡 
水 文 站 多 年 平均 径流 量 为 10.4x10* m ,最 大 值 出 现 
在 2017 年 ,为 15.92x10 m? ,最 小 值 出 现在 1997 年 ， 
H 5.13x10* m^, 1970—1990 年 呈现 出 明显 的 丰 枯 
年 交替 变化 ,1990 一 2005 年 基本 处 于 枯 水 年 份 ， 
2006 一 2020 年 基本 处 于 丰 水 年 份 。 从 线性 拟 合 情 
况 来 看 ,年 径流 量 总 体 呈 增加 趋势 ,变化 率 为 
0.0475x10* m .ar:。 从 Sa 滑动 平均 曲线 来 看 ,径流 
量 总 体 呈 上 上升- 下降- 上 升 的 波动 变化 趋势 。1990 
年 以 前 5 a 请 动 平 均 曲线 小 幅 波动 ,基本 位 于 平均 值 
上 下 ,表明 1970 一 1990 年 这 一 时 期 正义 峡 站 年 径流 
量 比较 稳定 且 基 本 接近 多 年 平均 值 ;1990 一 2007 年 
曲线 位 于 多 年 平均 值 下 方 , 说 明 这 段 时 间 内 正义 峡 
站 年 径流 量 较 少 ,其 中 ,1990 一 1994 年 曲线 下 降 趋势 
明显 ,1995 一 2003 年 较为 稳定 ;从 2004 年 开始 ,5 a 
滑动 平均 曲线 呈现 明显 的 上 升 趋势 ,表明 自 2004 年 


径流 量 /108ms 
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开始 年 径流 量 较 往年 增多 ,这 有 可 能 与 2000 年 开始 
实施 的 黑河 中 游 调 水 活动 有 关 。 

2.1.2 径流 突变 检验 ”正义 峡 水 文 站 年 径流 量 的 突 
变 检 验 如 图 3 所 示 。 图 3a 为 M-K 突变 检 验 , 图 中 UF 
统计 量 呈 现 先 上 升 . 后 下 降 .再 上 升 的 波动 趋势 , 且 
UF BE il t& fk 1970—1989 年 , 除 1980 4E. , 1986 年、 
1987 ^F. , 1988 年 以 外 ,其 余年 份 均 为 正 值 且 处 在 
0.05 显著 水 平 临界 线 (+1.96) 内 ,表明 正义 峡 水 文 站 
年 径流 量 在 此 段 时 间 呈 现 震 荡 上 升 趋势 ;1989 一 
2013 4E UF 统计 量 为 负 值 ,其 中 1997 一 2006 年 明显 
超过 0.05 显著 水 平 临界 线 (-1.96) ,表明 1989 一 2013 
年 正义 峡 水 文 站 年 径流 量 呈 现下 降 趋势 ,日 1997 一 
2006 年 下 降 趋势 显著 ;UF 和 UB 两 条 曲线 相交 于 
2014 年 左右 ,根据 M-K 突变 检验 规则 可 知 , 正 义 峡 
水 文 站 年 径流 量 在 2014 年 前 后 发 生 了 一 次 突变 。 
图 3b 为 滑动 :检验 ,图 中 统计 量 ; 在 2004 年 .2013 年 
左右 超出 给 定 P=0.05 显著 水 平 临界 线 , 其 中 ,最 小 
值 出 现在 2004 年 ,表明 正义 峡 水 文 站 年 径流 量 在 
2004 年 和 2013 年 前 后 分 别 发 生 了 显著 突变 ,其 中 
2004 年 突变 极 显 著 。 图 3c 为 Pettitt 检验 ,图 中 统计 
tit UT RÆ 1970 一 1990 年 变化 较为 平缓 ,1990 一 
2005 年 呈现 下 降 趋 势 ,2005 年 后 波动 上 升 , 并 超过 
P=0.05 显著 水 平 临界 线 。 根 据 Pettitt 突变 检验 中 最 
显著 的 可 能 突变 点 为 统计 量 UT 绝对 值 的 最 大 值 规 
则 ,图 3c 中 UT 绝对 值 的 最 大 值 出 现在 2005 年 且 超 
过 0.05 显著 水 平 临界 值 , 因 此 ,可 以 得 出 正义 峡 水 


一 - 正义 峡 径 流量 一 -5 a 滑动 --- 均值 …… 线性 (正义 峡 径 流量 ) 
y-0.0475x*9.1597 
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图 2 正义 峡 水 文 站 年 径流 变化 趋势 
Fig.2 Trend of annual runoff change at Zhengyi Gorge hydrological station 
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图 3 正义 峡 水 文 站 年 径流 突变 检验 


Fig.3 Test for annual runoff mutation at Zhengyi Gorge hydrological station 


文 站 年 径流 量 可 能 在 2005 年 前 后 发 生 了 一 次 显著 
突变 。 图 3d 为 累积 距 平 曲线 ,由 图 可 以 看 出 正义 峡 
水 文 站 年 径流 量 累 积 距 平 曲线 从 1970 一 1989 年 呈 
现 上 升 趋势 但 不 显著 ;从 1989 一 2004 年 整体 呈现 明 
显 的 单调 下 降 趋 势 ,上 且 下 降幅 度 较 大 ,说 明 在 此 段 
时 间 内 正义 峡 水 文 站 年 径流 量 较 少 ; 自 2004 年 后 ， 
累计 距 平 曲线 呈现 显著 的 单调 上 升 趋势 ,表明 年 径 
流量 增加 趋势 明显 ;年 径流 量 累 积 距 平 在 2004 年 前 
后 发 生 了 先 减 小 后 增 大 的 转折 变化 ,这 表明 正义 峡 
水 文 站 径流 量 可 能 在 2004 年 前 后 发 生 了 一 次 突变 。 

由 上 述 分 析 结 果 可 知 ,利用 M-K 突 变 检验 08 
动 ! 检 验 .Pettitt 检 验 和 累积 距 平 法 四 种 方法 对 正义 
峡 水 文 站 径流 量 的 突变 检验 结果 不 完全 一 致 ,但 其 


中 三 种 方法 均 表 明 突 变 点 在 2004 一 2005 年 前 后 。 
结合 图 2 中 年 径流 变化 趋势 综合 考虑 后 ,最终 确定 
正义 峡 水 文 站 径流 的 突变 时 间 点 为 2004 年 ,这 与 吴 
凯 等 1 的 研究 结果 一 致 。 

突变 前 后 径流 变化 情况 见 表 1。 由 表 1 可知, 突 
变 后 正义 峡 水 文 站 的 年 平均 径流 量 与 突变 前 相 比 
存在 较 大 幅度 的 增加 ,年 平均 径流 量 由 突变 前 的 
9.43x10* mm 增加 到 突变 后 的 12.51x10' m? ,增加 率 为 
32.7% ; 变异 系数 C, 由 0.28 2523 0.13 ,表明 突变 后 径 
流量 年 际 间 变异 程度 更 小 、 变 化 更 趋 稳定 。 
2.2 径流 变化 归 因 分 析 
22. 敏感 性 分 析 黑河 干流 中 游 区 域 面积 》 
2.56x10* km ,采用 正义 峡 水 文 站 多 年 实测 径流 量 资 


表 1 正义 峡 径 流 突变 前 后 对 比 
Tab.1 Comparison before and after the abrupt change in the runoff of Zhengyi Gorge 


K3 突变 点 突变 前 突变 后 
- 一 — - — - 变化 率 /% 
站 点 年 份 显著 性 径流 量 /10:m ”变异 系数 (C.) 径流 量 /10'm” ”变异 系数 (C.,) 

正义 峡 2004 0.05 9.43 0.28 12.51 0.13 32.7 


料 ,计算 径 流 深 民 ;利用 中 游 地 区 高 台 VE .张掖 三 
个 气象 站 的 实测 数据 ,通过 Penman-Monteith 公式 计 
算出 中 游 地 区 的 潜在 蒸发 量 EI7v。 

以 2004 年 作为 分 界 点 ,将 研究 期 (1970 一 2020 
年 ) 划 分 为 基准 期 (1970 一 2004 年 ) 和 变化 期 
(2005—2020 年 )。 根 据 多 年 平均 的 径流 深 (R) 、 降 
水 量 (P) 以 及 潜在 蒸发 量 (87,) 计 算 下 热 面 参数 nn、 
径流 系数 a(RAP) 和 干燥 指数 K(ETwP) ,再 通过 相关 
文献 中 的 公式 ”3 计算 得 到 不 同时 期 的 PETon 分 
别 所 对 应 的 弹性 系数 (se) ,计算 结果 见 表 2。 

由 表 2 可知, 正义 峡 水 文 站 变化 期 的 RP、ET。 
与 基准 期 相 比 均 有 所 上 升 , 上 升 速率 分 别 为 
38.02% .9.84% 和 10.01% ,径流 深 上 升幅 度 较 大 。 径 
流 系 数 a 较 基准 期 有 所 上 升 ,表明 变化 期 增加 的 降 
水 转化 为 径流 的 比例 比 基 准 期 高 ,干燥 指数 (K) 则 
基本 维持 不 变 ,符合 温带 干旱 区 的 气候 特征 。6, 为 
正 值 ,表示 降雨 量 (P) 与 径流 深 (R) 呈 正 相 关 关 系 ， 
当 P 每 增加 1 mm 时 ,R 的 增加 量 为 1.4 mm, He, JA AE 
准 期 的 1.49 减少 至 变化 期 的 1.40, 表 明正 相关 性 有 
减弱 趋势 ,这 可 能 与 潜在 蒸发 量 E7Tu 和 下 垫 面 参数 m 
与 径流 深 的 负 相 关 性 减弱 有 关 。 因 为 当 ET, 和 nn 分 
别 增 加 1 mm 或 1 时 ,径流 深 的 减少 量 由 基准 期 的 
0.49 mm 变 为 变化 期 的 0.40 mm ,或 由 1.79 mm 变 为 
1.57 mm。 说 明 径 流 变 化 不 单纯 受 某 一 变量 (P、 
ETon) 的 影响 ,而 是 多 变量 综合 影响 的 结 

弹性 系数 (e) 的 绝对 值 越 大 表明 径流 对 影响 因 
子 的 敏感 度 越 高 小 。 由 此 可 知 , 人 研究 期 内 正义 峡 水 
文 站 径流 对 下 热 面 参数 (n) 的 敏感 度 最 强 , 对 潜在 
蒸发 量 (ETv) 的 敏感 度 最 弱 。 
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2.2.2 归 因 分 析 通过 Budyko 水 热 耦合 平衡 方程 
计算 各 影响 因子 对 径流 变化 的 贡献 率 , 结 果 见 表 3。 

由 表 3 可 知 , 根 据 Budyko 模型 计算 得 到 的 径流 
深 变 化 量 (dR' ) 为 14.3 mm ,实际 径流 深 变 化 量 (dR) 
值 为 13.69 mm ,径流 深 变 化 量 的 相对 误差 为 4.46%， 
INF 5%" ,说 明 Budyko 模 型 适用 于 该 区 域 的 径流 
变化 研究 。 变 化 期 内 ,PET, 和 nn 对 径流 的 贡献 程度 
有 和 较 大 差异 ,气候 因素 (P 和 E7,) 对 正义 峡 水 文 站 径 
流 的 贡献 率 为 30.21% ,其 中 降水 引起 的 径流 深 的 变 
化 量 为 6.11 mm, 贡献 率 42.73% ,潜在 蒸发 引起 的 径 
流 深 的 变化 量 为 -1.79 mm ,贡献 率 -12.52% ;下 热 面 
引起 的 径流 变化 量 为 9.98 mm, 贡献 率 69.7996. £y 
上 可 知 , 正 义 峡 水 文 站 径流 量变 化 主要 归 因 于 下 热 
面 参数 的 变化 ,气候 因素 中 降水 量变 化 对 径流 量 的 
影响 相对 较 大 ,潜在 蒸发 量变 化 的 影响 次 之 。 


3 讨论 


3.1 气候 变化 对 径流 的 影响 

已 有 研究 表明 ,各 种 形式 的 大 气 降水 是 径流 的 
直接 补给 源 ,因此 ,降水 对 径流 产生 直接 影响 ;区 域 
气 侯 以 及 下 垫 面条 件 决 定 蒸 散发 程度 和 空间 格局 ， 
进而 对 区 域 径流 形成 以 及 径流 量 大 小 产生 影响 ; 近 
60 a 来 黑河 流域 径流 总 体 呈 现 出 增多 趋势 ,其 变化 主 
要 受气 温 与 降水 的 影响 ,其 中 降水 影响 更 大 ”i。 为 
了 进一步 分 析 黑 河中 游 正 义 峡 径流 量变 化 与 气候 
因子 之 间 的 关系 ,根据 已 知 资料 对 径流 量 与 降雨 
量 \ 气 温 以 及 潜在 燕 发 量 之 间 的 相关 性 进行 分 析 
(图 4)。 


表 2 正义 峡 不 同时 期 水 文 气象 变量 特征 


Tab.2 Characteristics of hydrometeorological variables in Zhengyi Gorge in different periods 


水 文 站 时 期 R/mm P/mm ET/mm n a K Ep &ET, £, 
正义 峡 ”基准 期 (1970 一 2004 年 ) 36.01 118.90 897.80 0.66 0.30 7.55 1.49 -0.49 -1.79 
变化 期 (2005 一 2020 年 ) 49.70 130.60 987.60 0.58 0.38 7.56 1.40 -0.40 -1.57 


性 系数 。 下 同 。 


表 3 正义 峡 径流 变 


注 :R 为 多 年 平均 径流 深 ;P 为 降水 量 ;ET 为 潜在 蒸发 量 ;n 为 下 垫 面 参数 ;a 为 径流 系数 ;KK 为 干燥 指数 ,sj、eBTh\es 分 别 为 PET,n 所 对 应 的 弹 


化 归 因 分 析 结 果 


Tab.3 Results of attribution analysis of runoff change in Zhengyi Gorge 


贡献 量 /mm 贡献 率 (C)/% 
水 文 站 基准 期 变化 期 dR'mm dRmm 6/mm 
P ET, P ET, n 
正义 峡 1970 一 2004 年 ”2005 一 2020 年 6.11 -1.79 9.98 14.30 13.69 0.61 42.73 =12.52 69.79 
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图 4 正义 峡 径 流量 与 气候 因子 的 相关 性 


Fig. 4 Correlation between runoff and climate factors in Zhengyi Gorge 


如 图 4a 所 示 ,1970 一 2020 年 黑河 中 游 正义 峡 径 
流量 变化 与 降水 量 P 变 化 趋势 基本 一 致 ,呈现 正 相 
关 关 系 , 但 相关 性 不 显著 。 从 基准 期 到 变化 期 径流 
量 随 降 水 量变 化 的 增加 率 从 0.0324 减少 到 0.0066 , 
说 明 变 化 期 降水 对 径流 的 贡献 在 减弱 ,影响 径流 年 
际 间 的 变化 可 能 还 与 其 他 因素 有 关 。 由 图 4b 可知， 
正义 峡 径 流量 与 气温 的 相关 性 呈现 微弱 正 相 关 , 从 
基准 期 到 变化 期 径流 量 随 气温 变化 的 增加 率 从 
1.0064 KEJ 0.8599 ,说 明 变 化 期 气温 对 径流 的 影响 
也 在 减弱 。 薰 散发 是 黑河 中 游 地 区 水 分 耗 散 的 主 


径流 量 愈 小 。 
结合 上 述 气候 因子 的 影响 分 析 ,本 文 表 3 计 算 
得 出 的 降水 量 (P) ,潜在 蒸发 量 (£7,) 对 正义 峡 径 流 
量变 化 的 影响 结果 ,与 相关 学 者 "的 研究 结果 一 
致 , 即 气 候 因 子 中 的 降水 对 黑河 中 游 正义 峡 径 流 演 
变 的 影响 占 主导 作用 。 
3.2 人 类 活动 对 径流 的 影响 

黑河 中 游 特殊 的 下 垫 面 和 地 质 构造 条 件 ,致使 区 
域 径流 的 形成 汇集 及 转换 可 能 还 会 更 复杂 一 些 ”。 
相关 研究 ”表明 ,下 热 面 参数 (n) 与 流域 的 地 形 、 


要 途径 之 一 ,受气 温 变 化 和 下 垫 面 等 因素 的 综合 影 ”土壤 性 质 .土地 利用 类 型 以 及 植被 覆盖 度 等 因素 有 


啊 ,夏季 气温 高 , 蒸 散 发 量 最 大 ,而 冬季 气温 低 , 蒸 
散发 量 最 小 ”。 在 一 定 的 气候 条 件 下 ,下 热 面条 件 
影响 蒸 散发 进而 影响 地 表 径 流 的 形成 及 其 年 内 分 
配 。 由 图 4c 可知 ,正义 峡 径 流量 与 潜在 蒸发 量 E7。 
的 相关 性 呈现 极 微弱 的 负 相 关 , 从 基准 期 到 变化 期 
径流 量 随 潜在 蒸发 量变 化 的 下 降 率 从 0.0026 减 小 
为 0.0024, 说 明 变 化 期 潜在 蒸发 量 对 径流 的 影响 也 
在 减弱 。 余 加 男 等 ”认为 降水 量 和 上 游 来 水 量 是 
黑河 中 游 主要 的 水 分 输入 变量 ,其 多 寡 决 定 了 中 游 
地 区 燕 散发 量 的 大 小 以 及 可 利用 水 分 的 多 少 。 在 
多 年 平均 的 时 间 尺 度 上 ,流域 水 量 输入 (降水 和 上 
游 来 水 ) 一 定 的 情况 下 ,气温 越 高 , 蒸 散发 量 愈 大 ， 


关 , 且 一 般 考 虑 到 在 较 短 时 间 内 流域 的 地 形 和 土壤 
性 质变 化 较 小 ,可 忽略 不 计 ,因此 ,流域 下 热 面 参数 
(n) 的 变化 主要 用 流域 的 植被 覆盖 度 及 土地 利用 类 
型 变化 来 反映 。 下 垫 面 条 件 变化 不 仅 影响 内 陆 水 
文 循环 和 流域 产 汇流 过 程 , 对 河川 径流 的 二 次 分 配 
过 程 也 产生 复杂 的 影响 。 为 进一步 分 析 人 类 活动 
对 河川 径流 变化 的 影响 ,本文 主要 采用 1981 一 2019 
年 黑河 中 游 地 区 的 归 一 化 植被 指数 (NDVI) 与 土地 
利用 面积 变化 指标 ,对 正义 峡 径 流量 年 际 变化 的 影 
响 进 行 分 析 , 如 图 5 所 示 。 

由 图 5a 可知 ,1981 一 2004 年 黑河 中 游 植被 ND- 
VI 呈现 下 降 趋 势 , 下 降 速 率 为 0.0039, 尤 其 在 
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Fig. 5 NDVI and land use area changes in the middle reaches of Heihe River 


1994—2004 ^F. , 随 着 草地 、 林 地 面积 的 减少 ,植被 
NDVI 呈 明显 下 降 趋势 。2004 年 之 后 草地 、 林 地 面 
积 有 小 幅度 增加 , NDVI 呈现 波 动 上 升 趋势 ,上 升 速 
率 较 小 , 仅 为 0.0023。 下 垫 面 植被 NDVI 的 这 种 变 
化 情况 与 正义 峡 径 流 突变 发 生 在 2004 年 相对 应 ,说 
明 进 入 变化 期 (2005 一 2020 年) 黑河 中 游 的 植被 覆 
盖 度 有 所 增加 ,也 表明 这 一 时 期 由 于 相关 水 资源 管 
理 政策 的 实施 ,水 资源 重复 利用 率 不 断 提高 ,加 之 
节 水 型 高 效 农业 的 推广 ,农业 用 水 量 减少 ,生态 植 
被 面积 有 所 增加 ,进而 对 正义 峡 径 流量 的 增加 产生 
正 影 响 。 但 从 图 $ 的 土地 利用 情况 可 以 看 出 ,变化 
期 耕地 面积 ,城乡 .工矿 .居民 用 地 面积 的 持续 增加 
将 会 影响 植被 生态 用 水 ,进而 影响 生态 植被 的 进 一 
步 恢 复 和 发 展 。 变 化 期 植被 NDVI 虽 呈 增 大 态势 ， 
但 未 达到 基准 期 的 植被 覆盖 水 平 ,表明 人 类 的 一 些 
活动 对 黑河 中 游 生态 植被 的 恢复 还 存在 一 定 的 负 
影响 。 

黑河 中 游 是 黑河 流域 水 土 资 源 开 发 利用 程度 
最 高 的 区 域 。 已 有 人 研究 ”显示 , 随 着 人 口 增加 和 
社会 经 济 发 展 ,人 研究 区 耕地 面积 城镇 化 建设 以 及 


水 利 设施 不 断 增加 , 极 大 地 改变 了 中 游 下 垫 面 条 
件 。 由 表 4 可 知 ,在 变化 期 内 (2005 一 2020 年 ) ,耕地 
面积 建设 用 地 面积 增长 率 分 别 为 10.33%、39.32%， 
AH GDP 的 增长 率 达 到 了 305.38%。 同 时 ,农业 用 
水 量 .工业 用 水 量 、 生 活用 水 量 则 呈现 下 降 趋 势 。 
这 与 张掖 市 进一步 深化 节 水 型 社会 建设 ,积极 落实 
《黑河 流域 近期 治理 规划 》 中 的 目标 任务 ,调整 产业 
结构 比例 ,统筹 地 表 水 和 地 下 水 资源 ,不 断 提 高 水 
资源 重复 利用 率 等 措施 有 关 。 黑 河中 游 也 属于 绿 
洲 灌溉 农业 区 ,农业 灌溉 用 水 约 占 中 游 总 用 水 量 的 
80% 以 上 ,在 变化 期 内 实际 耕地 灌溉 率 虽 增加 
5.1396 ,但 农业 用 水 量 降低 了 14.62% ,说 明 中 游 节 水 
型 高 效 农 业 的 发 展 成 效 明显 。 变 化 期 内 生态 用 水 
降幅 最 大 ,减少 了 60.31% ,说 明 中 游 地 区 为 了 保证 
农业 灌溉 用 水 ,挤占 了 区 域 生 态 用 水 ,这 也 可 能 是 
变化 期 内 中 游 生态 植被 覆盖 度 (NDVI) 增 速 缓慢 的 
原因 之 一 。 黑 河中 游 地 势 平坦 ,河道 落差 小 ,在 185 
km 的 天 然 河 道上 有 许多 不 同 规模 的 蕃 水 引水 工程 ， 
但 无 水 电站 工程 ,这 些 攻 水 引水 工程 在 为 工农 业 生 
产 及 居民 生活 提供 水 源 的 同时 ,也 不 同 程度 改变 了 


表 4 黑河 中 游人 类 活动 指标 的 年 际 变化 


Tab. 4 Interannual variation of human activity indicators in the middle reaches of Heihe River 


fy 耕地 面积 / 实 Mk 地 人 均 GDP 建设 用 地 农业 用 水 ”工业 用 水 生活 用 水 ”生态 用 水 总 用 水 量 
km? 灌溉 率 /96 /10: 元 1 积 /km /10° m° /10* n /10: m° /A0* m° /A0* m° 
2000 1849.67 72.4 0.32 176.73 16.64 0.51 0.37 1.42 18.94 
2005 2003.50 85.96 0.93 189.32 15.32 0.66 0.48 1.31 17.77 
2010 2133.49 84.65 1.77 200.62 15.67 0.30 0.40 1.22 17.59 
2015 2178.12 83.91 3.06 224.57 15.76 0.22 0.41 0.78 17.18 
2020 2210.45 81.55 3.77 263.77 13.08 0.10 0.36 0.52 14.06 
变化 期 (2005 一 10.33 -5.13 305.38 39.32 -14.62 -84.85 -25.00 -60.31 -20.88 


2020) 变 化 率 /% 
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中 游 原 河道 的 径流 特征 和 地 表 径 流 分 配 过 程 ,间接 ” 感性 最 弱 。 


影响 区 域内 土地 利用 和 生态 植被 的 空间 格局 。 

另外 ,为 了 遏制 黑河 下 游 生态 退 化 的 趋势 , 自 
2000 年 8 月 黑河 流域 实施 了 《黑河 干流 水 量 分 配方 
案 》, 黑 河中 游 每 年 在 营 落 峡 多 年 平均 来 水 15.8x10’ m 
的 情况 下 分 配给 正义 峡 9.5x10 m? B FKE". E 
至 2020 年 ,累计 进入 正义 峡 断 面 水 量 逾 225x10* m^, 
占 到 营 落 峡 断 面 来 水 总 量 的 60.5%, 极 大 地 改善 了 
黑河 干流 下 游 的 生态 环境 2 。 在 干旱 地 区 , 基 流 
是 河川 径流 的 一 部 分 , 基 流 在 维持 河道 水 源 稳定 性 
和 脆弱 荒漠 生态 系统 中 发 挥 着 重要 作用 。 廉 次 康 
等 (3 分 析 发 现 正义 峡 水 文 站 1957 一 2017 年 期 间 年 
基 流 指数 值 达到 了 0.7571, 调 水 前 后 正义 峡 基 流 指 
数 变化 与 地 表 径 流 变 化 不 尽 同 步 ; 调 水 前 , 基 流 呈 
现 减 小 趋势 ,地 表 径流 呈现 较 大 的 减 小 趋势 ; 调 水 
后 , 基 流 呈现 增加 趋势 ,而 地 表 径 流 呈 现 更 大 的 增 
加 趋势 。 可 见 ,黑河 中游 调 水 活动 ,在 增加 正义 峡 
径流 量 的 同时 ,无 疑 也 会 对 中 游 河 道 的 径流 分 配 、 
利用 以 及 区 域 生 态 环境 变化 产生 较 大 的 影响 。 

综 上 所 述 ,正义 峡 径 流量 的 年 际 变 化 与 中 游 地 
区 的 气候 条 件 及 人 类 活动 息息相关 。 在 自然 条 件 
与 人 类 活动 作用 下 ,黑河 中 游 地 表 水 与 地 下 水 相互 
转化 频繁 ,降水 .蒸发 .河道 渗 漏 灌溉 引 排 水 .地 下 
水 溢出 ,形成 了 干旱 区 特殊 的 水 资源 演变 模式 。 加 
强 对 中 游 正 义 峡 径流 演变 规律 的 研究 ,深化 气候 变 
化 和 人 类 活动 对 径流 影响 程度 的 认识 ,对 于 黑河 中 
游 水 资源 优化 配置 和 管理 有 具有 指导 意义 。 


4 结论 

本 文 以 正义 峡 水 文 站 实测 年 径流 量 数据 为 依 
据 ,对 正义 峡 水 文 站 1970—2020 年 径流 变化 及 归 
进行 了 分 析 , 得 出 以 下 主要 结论 : 

(1) 1970 一 2020 年 正义 峡 径 流量 年 际 间 波动 变 
化 ,但 总 体 呈 增加 趋势 ;正义 峡 径 流量 在 2004 年 发 
生 突 变 , 变 异 系数 C, 由 突变 前 的 0.28 变 为 突变 后 的 
0.13 ,表明 发 生 突变 后 正义 峡 年 径流 量变 化 呈 弱 变 
异 ,径流 变化 更 趋 稳定 。 

(2) 变化 期 内 (2005 一 2020 年 ) 正 义 峡 径流 量 对 
降水 .潜在 蒸发 和 下 垫 面 参数 的 弹性 系数 分 别 为 
1.40 -0.40 -1.57, 表 明正 义 峡 径流 量 对 中 游 下 垫 面 
变化 最 为 敏感 ,降水 变化 次 之 ,潜在 蒸发 变化 的 敏 


(3) 黑河 中 游 的 降水 、 潜 在 蒸发 和 下 垫 面 参数 
对 正义 峡 径流 的 贡献 率 分 别 为 42.739% -12.52% 、 
69.79% ,表明 人 类 活动 引起 的 黑河 中 游 下 垫 面条 件 
变化 是 影响 正义 峡 径流 量 演变 的 主要 原因 。 但 正 
义 峡 径流 演变 不 单纯 是 某 一 变量 (PET、m) 的 影 
响 , 而 是 多 种 变量 综合 影响 的 结 
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Runoff trend and attribution analysis of the Zhengyi Gorge in the middle 
reaches of the Heihe River 
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(1. School of Environmental and Municipal Engineering, Lanzhou Jiaotong University, Lanzhou 730070, Gansu, 


China; 2. Key Laboratory of Yellow River Water Environment in Gansu Province, Lanzhou 730070, Gansu, China) 


Abstract: This study utilized measured runoff data from the Zhengyi Gorge hydrological station in the middle 
reaches of the Heihe River to analyze the trend and mutation years of the runoff sequence from 1970 to 2020. 
Various methods, including linear regression, the Mann- Kendall mutation test, the sliding 7' mutation test, the 
Pettitt test, and cumulative offset verification, were employed for the analysis. Additionally, an attribution 
analysis of the runoff changes at the Zhengyi Gorge hydrological station was conducted using the Budyko water- 
energy coupled balance equation. The results indicate the following: (1) Over the study period, the runoff at 
Zhengyi Gorge showed fluctuations, alternating between periods of abundance and contraction but showed an 
overall increasing trend. The runoff experienced a mutation in 2004, resulting in a 3.08 x 10° m? increase in 
average annual runoff volume, representing a growth rate of 32.7%. (2) In the period after the mutation (2005- 
2020), the elastic coefficients of runoff in response to precipitation, potential evaporation, and underlying surface 
parameters were 1.40, —0.40, and — 1.57, respectively. The contribution rates of each factor to runoff were 42.7395, 
—12.52%, and 69.7995, respectively, indicating that runoff is most sensitive to underlying surface changes and 
that precipitation has a more significant impact on runoff than potential evaporation among climate factors. (3) 
Under certain regional climatic conditions, the changes in the underlying surface of the middle reaches caused by 
human activities, such as vegetation cover, land use, and watershed water transfer, are the main reasons for the 
changes in runoff at Zhengyi Gorge. The research findings can provide a scientific basis for river basin 
management departments to formulate water resource allocation and utilization plans. 


Keywords: runoff change; attribution analysis; coupled energy-water balance equation; climate change; human 
activity; the middle reaches of Heihe River 


